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In these notes, the application of Feynman's sum-over-paths approach to 
thermal phase transitions is discussed. The paradigm of such a spacetime 
approach to critical phenomena is provided by the high-temperature ex- 
pansion of spin models. This expansion, known as the hopping expansion 
in the context of lattice field theory, yields a geometric description of the 
phase transition in these models, with the thermal critical exponents be- 
ing determined by the fractal structure of the high-temperature graphs. 
The graphs percolate at the thermal critical point and can be studied 
using purely geometrical observables known from percolation theory. Be- 
sides the phase transition in spin models and in the closely related </> 4 
theory, other transitions discussed from this perspective include Bose- 
Einstein condensation, and the transitions in the Higgs model and the 
pure U(l) gauge theory. 

1. Introduction 

Feynman's spacetime approach to quantum mechanics provides a mar- 
velously intuitive, yet computationally powerful entry to the quantum 
worldMJ In this sum-over-paths approach, transition amplitudes are com- 
puted by summing over all possible trajectories a particle can take. 
Although it is most extensively applied to nonrclativistic quantum 
mechanics P for which it was originally designed, the worldline app roach 



also found applications in the context of relativistic systems!^^ (For 
overviews and references to a host of applications in high energy physics, 
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see Refs. [5J [51 [S].) In these notes, we explore its use in the theory of 
critical phenomena. 

The possibility to adapt Feynman's sum-over-paths approach to describe 
phase transitions derives from particle-field duality. Featuring the particle 
content of the system, the spacetime approach provides a completely equiv- 
alent alternative to the quantum field theoretic description. Although by 
no means the only successf ul metho d,^ field theory is frequently used to 
describe critical phenomena) ^ ^ 1 1 1 J12I whether the fluctuations are thermal 
as in phase transitions taking place at finite temperature, or quantum as 
in zero-temperature transitions. Particle-field duality then opens the pos- 
sibility for a spacetime approach to critical phenomena. In the context of 
classical critical phenomena, the line objects described by the field theory 
are current lines. 

Th ough not formulated as such, the high-temperature expansion of spin 
models^lfil probably provides the first example of the spacetime approach 
to critical phenomena (see below). A particular intuitive example is offered 
by the nonlinear O(N) spin model defined on a lattice in the limit N — > 0, 
where the high-temperature graphs form random self-avoiding walksff^l I n 
the early 1950ies, FeynmarP^ applied his approach to the superfluid phase 
transition in liquid 4 He, and showed that, in this picture, Bose-Einstein con- 
densation arises when the worldlines of individual particl es combine to fo rm 
large lines threading the entire system. In the late 1970ies p^^^^^^^^^ lthe 
approach was applied to phase transitions in relativistic lattice models pos- 
sessing line-like topological defects, thus mapping these models onto the 
statistical mechanics problem of random loop configurations. The line de- 
fects were pictured as forming a grand canonical ensemble of fluctuation 
loops of arbitrary shape and length-a so-called loop gas. In the 1980ies, 
the statistical properties of such loop gases were further investigated in the 
context of cosmic strings and their role in phase transition s in th e early 
universe, with the loop distribution given a central positioning as we n 
as in the context of deconfinement phase transit ions, w here percolation ob- 
servables were introduced to characterize themPSUl More recently, these 
ideas were pursued further with the goal to arrive at a complete quantitative 
understanding of phase transitions i n terms of the geometr ical properties 
of the loop gas under 

considerationP26|27|28|29|30|31|32| 

The resulting 

geometric description of critical phenomena ties together various strains 
of theoretical physics: statistical mechanics, quantum field theory, polymer 
physics, and percolation theory-all to be discussed below. 

The plan of these notes is as follows. In Sec. El the particle content of 
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the linear 0(A) symmetric quantum field theory, defined on a hypercu- 
bic lattice, is studied in the spirit of Feynman's sum-over-paths approach, 
starting from the noninteracting theory in Sees. 12.11 and 12.21 Interactions 
are perturbatively included in Sec. 12. 31 where also the spacetime interpre- 
tation of Feynman diagrams is pursued. In Sec. 12.41 the strong-coupling 
limit, where the field theory reduces to the nonlinear 0(A) spin model, 
is considered. The high-temperature representation of the nonlinear model 
is shown to provide a purely geometric representation, directly connected 
to Feynman's spacetime approach for the linear model. As a case study, 
the high-temperature representation of the Ising model, corresponding to 
A — 1, on a square lattice is detailed in Sec. 12.51 

In Sec. the critical properties of the 0(A) universality class are ex- 
amined from the perspective of particle trajectories or high-temperature 
graphs. In Sec. 13.11 their fractal structure is investigated, while in Sec. 13. 2\ 
the proliferation of these lines is shown to drive the 0(A) phase transi- 
tion. In Sec. 13.31 the fractal structure of the worldlines or high-temperature 
graphs is related to the critical behavior. The general relations are first ap- 
plied to self-avoiding random walks, corresponding to the N — ► limit, in 
two dimensions (Sec. I3.4j ). and to arbitrary —2 < A < 2 fSec. 13. 5|) also in 
two dimensions, where many exact results are available for comparison. 

In Sec. 21 the high-temperature representation of the Ising model is 
investigated by means of Monte Carlo simulations. The Metropolis update 
used to generate the high-temperature graphs is introduced in Sec. 14. 1[ 
while the numerical results are summarized in Sec. 14.21 

In Sec. 03 the conclusions for the 0(A) universality class are extended to 
charged systems (Sec. 15. ip and to Bose-Einstein condensation (Sec. 15. 2p to 
demonstrate the generality of the spacetime approach to phase transitions. 

Finally, in Sec. [SI the worldlines or high-temperature graphs are reinter- 
preted as line defects in two (Sec. 16. ip , three (Sec. l6.2p . and four (Sec. 16. 3p 
dimensions, respectively, and the so-called dual field theories, which feature 
these line-like configurations as topological defects, briefly discussed. 

2. Lattice Field Theory 

Central to our discussion is the 0(A) symmetric relativistic quantum field 
theory formulated on a hypercubic lattice in d Euclidean spacetime dimen- 
sions, with the time coordinate analytically continued to the imaginary 
axis t. As an illustration, the cubic lattice in Fig. Q] represents discretized 
spacetime in one time and two space dimensions x\ and X2- The theory, 
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Fig. 1. Cubic lattice representing three-dimensional discretized spacetime. 



describing self- interacting scalar particles of mass m, is specified by the 
(Euclidean) partition function 

-S 



Z = Tre 



where S is the lattice action 



5* = 



,d \ " 



1 v - 2 m2 2 9 

^2 2^ W»*+«»A - Vx) + ~V X + ^ 



(1) 



(2) 



with g the coupling constant and a the lattice spacing. The real scalar 
field ip x defined on each lattice site x of the spacetime box has N com- 
ponents ip x — ip% = (fl, ip x , ■ ■ ■ ,<Px), with the index a = 1, 2 • • • , N la- 
beling the field components and Lp\ = (ip x ■ <p x ) 2 , where the dot product 
implies a summation over the field components ip x ■ ip x = 'Y^=\ l P x '<P x - 
Lattice coordinates, representing discretized spacetime, are specified by 
x — — (xi, X2, • • • , Xd), with fi = 1, 2, • • • , d and Xd = t denoting the 
imaginary time coordinate, while fi denotes the unit vector pointing in the 
/^-direction. In the continuum limit a — > 0, the lattice action reduces to 



S 



d«x ] - + — + 



(3) 



where f(x) stands for the field defined in continuous spacetime. Unless 
otherwise indicated, natural units h = c = 1 are adopted throughout. The 
trace Tr in Eq. ([1]) stands for the sum or integral over all possible field 
configurations: 



(4) 



In the continuum limit, where the lattice spacing is taken to zero a — » 0, 
the right hand defines the functional measure denoted by J Dip. 
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For numerical simulations, a more convenient form of the lattice action is 

obtained by casting Eq. ((2]) in terms of dimensionless fields and parameters 

fool 

defined via^ 

a d -\l = 2K4> 2 n (5) 

- A ~ d 9 = 6^ (6) 



2 2 

m a 



K 2 
1 - 2\N 



2d. (7) 



K 

The action then takes the form of an O(N) spin model 

S = -K <A» • <t>n< + E + A E « -Nf. (8) 

(n,n') n n 

Each site on the hypercubic spacetime lattice is now specified by the vector 
n = rifx = (ni,ri2,-- - , rid) with integer components — x^/a and the 
sum X){n ri) extends over all nearest neighbor pairs. In terms of these new 
dimensionless variables, the action is independent of the lattice spacing a. 
The partition function Z now reads 

Z= D M (0)exp \K ^ (9) 

J \ (n,n') J 

with the on-site measure 

/OM0)= /n^e-^-^-") 2 . (10) 

In the limit A — > oo, the field theory reduces to the standard O(N) spin 
model, with a "spin" variable <j> n of fixed length, § 2 n = N, located at each 
site of the spacetime lattice. 



2.1. Noninteracting Theory 

We first consider the limit A — > in Eq. ([5]) , corresponding to the noninter- 
acting field theory. The rescaling factor K introduced in Eqs. ©-([7]) then 
becomes 

and the action reduces to the quadratic form 

So = ^ 4> n ■ A n<n /(f> n i, (12) 
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with 



^n,n' — $n,n' — K } j & n ,n'+\i (13) 

where the sum JZ±g extends over the positive as well as the negative di- 
rections. In matrix notation 



A = I - K H, 



(14) 



where / is the identity matrix and H denotes the so-called hopping matrix 
whose elements H n ^ n i are unity if the two lattice sites n and n' are near- 
est neighbors, and zero otherwise. Physically, H describes the hopping, or 
propagation of a particle from one lattice site, representing a spacetime cell, 
to an adjacent one, with each hop carrying a weight K. Since the theory 
is noninteracting, a particle is free to hop at random to any of its nearest 
neighbors without restriction. It can and will in general revisit sites pre- 
viously visited and can even hop onto a site already occupied by another 
particle. The particle trajectories are therefore phantom worldlines that can 
freely intersect and share bonds of the lattice without energy penalty. 
The inverse of A determines the lattice correlation function: 



as can be read off from the action (fT2|l . From the definition of the hopping 
matrix H it follows that H n ^ n i raised to the power b gives the number of 
paths Zb(n, n') of b steps starting at site n and ending at site n' . In this so- 
called hopping expansion^^ the correlation function is therefore obtained 
by summing over a ll po ssible worldlines of arbitrary many steps joining the 
endpoints n and n'!^ 



Each path carries a weight K b according to the number of steps b taken. 
This geometric representation of the correlation function as a sum over 
worldlines constitutes Feynman's spacetime approach^ to fluctuating fields 
on the lattice. The equivalence of the two approaches is known as particle- 
field duality. 

Because of translational invariance, the number of particle trajectories 
Zb{n,n') depends only on the relative coordinate Zb(n, n') — zt(n — n'). It 




(15) 




(16) 
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satisfies the recurrence relation 



z b+ i{n) = ^2 z b(n + A)> 



(17) 



with the initial condition Zo(n) — 5o,n- This relation is easily solved by 
going over to Fourier space 



z b (n - ri) = 



K r\ d n 

7T (^r 



with q ■ n = q^n^ = q\n\ 
variable, yielding 



Zb(q), (18) 
qdjid and g M the dimensionless momentum 

Zb+i(q) = cos( qil )z b (q). (19) 
With z (q) = 1 as initial condition, this gives for z b (n — n') 



2^cos(<? At ) 



(20) 



In the limit where the lattice spacing a is small but still finite, so that the 
particles still take discrete steps, we have 



2K 22 cos (<?/j 



1 - q 2 /2d 



1 + m 2 a 2 /2d 
-> exp [-(q 2 + m 2 a 2 )b/2d] 

The lattice correlation function (JTHJ) then becomes 

1 ~ d dg 



(21) 



3a, p 



2 J n (2v) d 



d/2 



.q-(n-n') e -(q 2 +m 2 a 2 )b/2d 



S «A E(^) 

6=1 v 7 



exp 



d (n - n') 2 
2 6 



(22) 



showing that the end-to-end vector n — n' is Gaussian distributed with an 
extra Boltzmann weight exp(— m 2 a 2 /2d) associated with each step taken. 
Because of this bond fugacity, long worldlines are exponentially suppressed 
as long as the mass term in the action is positive, m 2 > 0. Notice that in 
the last line of Eq. (|22|) . the lower bound on the summation is replaced by 
6=1. This is justified because the 6 = term always vanishes for n^n'. 
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A more familiar expression for the correlation function is obtained from 
Eq. (fT6| by first summing over the length b of the particle trajectories. With 
Zb(n,n') given in Eq. ([20| . this leads to 

1 r d d q e iq < n - n '^ 



2j_, (2^1-2/^cos^) 

0a ^2Kj_, (27r)*2£ M [l-cos(g M )]+m 2 a 2 - 1 Jj 

For the original fields i^, which are connected to the rescaled fields <fi n 
through Eq. |(SJ), this translates into the standard form of the correlation 
function of a free field theory on a hypercubic lattice, 

/ g\ d d k e ik-{x-x') 

= *aj* ^ 2 ^ {l _ cosM]+m ^ (24) 

where the lattice sites are now labeled again by = n^a, and k^ = q^/a 
is the momentum variable. 

In the limit a — > 0, the correlation function (pM)) reduces to the standard 
continuum expression 

/ a \ f°° d d k e lk -( x ~ x> ) 

The momentum integral can be easily carried out by introducing the 
Schwinger proper time representation^ of the integrand, which is based 
on Euler's form 

1 1 / "° d Ve-<S (26) 



A* T(z) J s 
where T(z) is the Euler gamma function. Specifically, 

d d k 



= <W/ q ds J A± e ^) e -s(k^) 



d/2-1 

i in i 

' (27r) d / 2 

with K d / 2 -i a modified Bessel function. The integral in the second line is 



S ^7^m(u-Zn) (27) 



the continuum limit of the lattice expression (|22p , with the sum over the 
number of steps b taken by the particle replaced by an integral over the 
Schwinger proper time parameter s, with the identification 

ba 2 /2d=s. (28) 
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The parameter s is related to the proper time r p of the trajectory traced 
out by the point particle, for which 

2 

= 1- (29) 



d^(rp) 



dr n 



To make this evident, we write the integrand in the second line of 
Eq. ([27|) as a Feynman path integraff^ 



Airs 



d/2 



e -(x-x') 2 /is g-sm 2 



x(s) = 



Dz(s') exp ^ - / ds' [±±V) + m 2 ] 
where the so-called Wiener measure is defined by 

x(s)= x' b—1 »oc 

Dx(s') = lim A b TT / d Xl , 



lx(0)=x 

with the normalization factor 



A = 



2na 2 



d/2 



(30) 



(31) 



(32) 



Here, it is understood that the limit b — > oo is taken concurrently with the 
continuum limit a — > in a way such that ba 2 — > const. The classical action 
of the point particle in Eq. (|30[) . tracing out its trajectory in spacetime is 
the continuum limit of the lattice action present in Eq. (|22[) : 



w. = IE 

i=l 
6 2 

a 



2d 



1 / Xi—i X\ 

4 V a 2 /2d 



2d 



ds' [\x 2 {s')+m 2 



(33) 



where in the last step the identification (f2"5| is used. Written as a path 
integral, the correlation function (f27l) thus reads 

/>oo /*x{s)—x' 

<pM)=S a ,p / ds / Dz( S ')e- w °. (34) 

' JO Jx(0)=x 

It forms the continuum limit of the sum over all possible particle trajectories 
on the spacetime lattice in Eq. (fT6|) . Note that the proper time parameter 
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s can take any positive value, in agreement with the lattice representation 
(| 16(1 of the correlation function where the worldlines can be arbitrarily long. 
The Boltzmann factor exp(— sm 2 ) exponentially suppresses, however, large 
proper times as long as m 2 > 0. When we set 

s = T v /2m, (35) 

the continuum action (|33[) yields m times the arc length of the path, 

ds' [j± 2 (s') + m 2 ] = mr p , (36) 

o 

which is the standard form of the classical action of a relativistic free scalar 
particle of mass m. This demonstrates that s can indeed be thought of as 
the proper time. 

For later use, we record the second line in Eq. (|27|) with Planck's con- 
stant and the speed of light reinstalled: 

POO / 1 \ 

(tf^) = <W C j( d S (-L- j e -(cW )/4 Ks (37) 

while the continuum action © becomes 



S = / drd^x ip ( -h 2 - h 2 V 2 + m 2 c 2 ) if, C-!N) 



where now the spacetime coordinates are denoted by x M = (x, cr), with r 
the imaginary time and x = {x\, • • • , Xd-i) the spatial coordinates. 

Besides the number of worldlines Zb(n, n') on the spacetime lattice con- 
necting the lattice sites n and n' in b steps, also the probability Pb(n,n') 
for a particle to move from n to n' in b steps will become important below. 
The probability is given by the ratio 

P b (n, n') = z b (n, n')/z b (39) 

of Zb(n, n') and the number of paths of b steps starting at n and ending at 
an arbitrary lattice site, 

z b = z b (n, n). (40) 



The denominator in Eq. ([39]) is included so that Pb is a genuine probability, 
taking values between and 1. To evaluate it, imagine the particle taking a 
first step to one of its nearest neighbors n±/2. It then has only n— 1 steps still 
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available to reach the final destination n' with probability P b -i(n ± p,, n'). 
That is, P b satisfies the recurrence relation: 



Pb(n,n') = — Y^Pb-i{n + jl,n') 



(41) 



where the factor 2d denotes the number of nearest neighbors on a d- 
dimensional square lattice. When Pf,_i(n, n') is subtracted from both sides 
of Eq. (|4ip. the difference equation 



P b (n,n')- P b -x(n,n') = ^ ^ [P b -i(n + £, n') - P b -i (n,n')] (42) 



follows, with the boundary condition P ( n 7 n ') = ^n,n'- ^ s before, the con- 
tinuum limit is taken by simultaneously letting a — > and & — > oo, such 
that ba 2 — > const. The difference equation then turns into the differential 
equation 



d b P b (n,n') = —V 2 P b (n,n') 



with the celebrated solution 



d/2 



d(n — n') 2 
2 b 



(43) 



(44) 



For the probability density p b (x,x') = P b (n,n')/a d this gives the Gaussian 
distribution 



Pb(x,x') 



2nba 2 



d/2 



exp 



d(x-x') 2 



2 ba 2 



(45) 



Each dimension the particle is free to roam contributes a factor \J d/2nba 2 
to the prefactor. 



2.2. Loop Gas 

A similar geometric representation in terms of particle trajectories as found 
for the correlation function can be given for the partition function. Consider 
the partition function of the free theory: 



*o=n/< 



-So 



(46) 
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with the lattice action Sq given by Eq. (fT!?|) . This integral generalizes the 
standard Gaussian integral 

n \ 

= (2^)"/ 2 det- 1/2 (M), (47) 



[ jjda^exp J -~ ^2 XiMijXj J = (27r 
J »=i \ ij=l / 



with M a symmetric positive-definite n x n matrix, and det(M) its deter- 
minant. It can therefore be evaluated exactly, with the result 

lnZ = lnDet _JV/2 (A) = — TrlnA, (48) 

where Tr denotes the functional trace. In deriving this, irrelevant prefactors 
are ignored and use is made of the identity 

lnDct(A) = TrlnA, (49) 

relating the determinant and trace of an Hermitian operator A. The relation 
can be easily checked explicitly for a n x n matrix. With A = I — KH as 
before, In Zq can be expanded in the hopping parameter K as 

N °° 1 

lnZ = -L d Y,TZb(a)K b (50) 

b=l 

where Zb(a) denotes the number of worldlines of b steps, starting at a given 
lattice site n and returning to one of its nearest neighbor sites n ± fi after b 
steps (see Fig. [2]). The lattice spacing a serves here as a microscopic cutoff, 
so that Zb(a) rather than Zb (0) appears when closing open paths. Since 
a worldline can start at any site on the spacetime lattice, the right hand 
contains a factor L d , denoting the total number of lattice sites, and lnZ is 
extensive. The factor 1/26 appearing in Eq. ([50| reflects the possibility to 
trace out a closed path starting from any site along the path and going in 
either direction round it. The factor N reflects the degeneracy associated 
with a closed worldline. Written as a sum over single loops of arbitrary 
length and shape, Eq. (|50|) reads 

lnZ = N K b , (51) 

single 
loops 

where the loops are considered to have no orientation, whence the absence 
of the factor |, and no longer rooted, whence the absence of the factor L d . 

For a noninteracting theory, Zb(a) can be obtained simply from Eq. 120|) 
by putting the end-to-end vector n — n' to zero there. Equation ([50l then 
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Fig. 2. A Brownian random walk on a square lattice returning to a site adjacent to 
its starting position n. As indicated, the particle took 27 steps in total. Observe that 
various bonds are traversed more than once. 



reads explicitly 



—L d r 

2 J-„ (2?r) d 



6=1 



Repeating the same steps leading to Eq. ([22]) . one finds 

r * d d q 1 



2 ^J_ n (2n) d b 

■ l(J_\ d/2 



)6/2rf 



2 ^ 6 V 2tt6 / 

6=1 v 7 



(52) 



(53) 



It follows that the logarithm of the partition function has the 

form™ 



with the loop distribution 



It, ~ 6~ r e- efc . 



(55) 



denoting the average number (per lattice site) of times a closed worldline of b 
steps occurs on the lattice. The parameter 6 oc m 2 physically determines the 
line tension of the loops, while the algebraic factor in the loop distribution, 
characterized by the exponent r = d/2 + 1, is an entropy factor, giving a 
measure of the number of ways a loop of b steps can be embedded in the 
lattice. 
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A more familiar expression for In Z Q is obtained by again first carrying 
out the sum over the loop length b in Eq. (|52|) : 



2 J_ n (2*)* 



l-2if^cos(g M ) 



s ln 2^[l-cosMl 



" n / 77; — rr hJ 2 V [1 - cos (fc M a)] + mV i , (56) 



2 J-,/a (2tt) 

where in the last step an irrelevant additive factor is ignored and f2 = (La) d 
is the volume of the spacetime box. This is the standard representation 
of the logarithm of the partition function describing a free field theory 
on a hypercubic lattice. In the continuum limit a — > 0, the momentum 
integral can again best be carried out by using the Schwinger proper time 
representation this time of the logarithm, 

In(a)= Urn- - / — e" sa , (57) 

where the first term is an irrelevant diverging contribution. Apart from 
further irrelevant additive constants this leads to 

(-81 + m 2 ) 

0^\d V 



= In Det 


N/2 (- 




/•CO 


-£\ 






Loo ( 


= ^nf 


00 ds 






2 Jo 


s 



m 2 ) 



1 x d/2 

I 1 

47TS / 



(58) 



where the operator — d 2 
Eq. (flU]) . The right hand of Eq. (|55|) forms the continuum limit of the lattice 
expression in the last line of Eq. (|53[) with the number of steps b taken by 
the particle assuming the role of the Schwinger proper time parameter s: 
ba 2 /2d= s as in Eq. ([28]) . 

On account of the path integral (|30|) with x = x', the last two factors in 
the integrand at the right hand of Eq. ([58]) can be written as a path integral 
over closed worldlines 



1 = j>Vx(s') ex p{-^ ds ' [ji 2 (s')+m 2 ] 



(59) 



where the notation § T)x(s') refers to closed paths. In this way, IuZq ex- 
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pressed as a path integral readd^ 



lnZ = lnDet _JV/2 (-dir. + /// 



2^ 

■ ■;- -in 



— n r °° ds 

~2 







j>T>x(s r ) exp ^- J ds' [±± 2 (s') + to 2 ] j , (60) 

where the argument of the exponential function is given by (minus) the 
classical action Q33[) of the point particle tracing out its trajectories in 
spacetime. An easy mnemonic for remembering the main ingredients of this 
formula is to interpret — d 2 + to 2 as the Hamilton operator describing the 
motion of a nonrelativistic particle of mass M = h in a constant potential 
V = to 2 in d dimensions. With the Schwinger proper time s interpreted as 
Newtonian time, Wq is the corresponding nonrelativistic action, where it is 
recalled that in natural units h=l. 

The first line in Eq. (|58|) is sometimes represented by the one-loop Feyn- 
man diagram 

lnZ = O- (61) 

In the spacetime approach, this diagram stands for all the terms in the sum 
(|50|) . each consisting of a single worldline starting and ending at a given site 

n. That is, the Feynman diagram at the same time denotes the topology of 

1041 

the particle trajectories contributing to lnZn,. 

The partition function, obtained by exponentiating Eq. (|50|) . can be 
written in terms of Feynman diagrams as 

Zo = i + + ^00 + ^000 + -- - (62) 

Each diagram involves a separate integration over loop momentum k^. With 
the right hand of Eq. (|62p taken as picturing the topology of the particle tra- 
jectories contributing to the partition function, it follows that Z describes 
a loop gas-a grand canonical ensemble of fluctuation loops of arbitrary 
shape and lengthPSl Being a noninteracting theory, the closed worldlines 
can freely intersect themselves or other worldlines and even share bonds. 
On the spacetime lattice, the loop gas can be expressed as a sum over- 
all possible tangles of closed paths 

Z Q = ]T K b N l , (63) 

loops 

with b and I denoting the total number of occupied bonds and separate 
loops forming a given tangle. In Eq. (|63|) . the loops are no longer consid- 
ered to be rooted, i.e., they can start and end at any spacetime coordinate 
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in the system. In comparison with rooted closed paths, i.e., closed paths 
all starting at a given spacetime coordinate, this yields a spacetime vol- 
ume factor f2. In the continuum, Eq. (|63|) translates into the path integral 
representation 

" e- w °\ (64) 



as follows from Eq. ([60)) . Here, Wn denotes the extension of the single 



particle action ([33]) to I particles: 



< )= E [^«)+m 2 ]. (65) 
1 Jo 



2.3. |0| 4 Field Theory 

We next perturbatively include the interaction in the theory @ . At the level 
of the partition function l[T]). this boils down to expanding the interaction 
term of the action in a Taylor series 



^(|) 2 « 2d E^ + -- (^) 



The ter ms in this expansion can all be represented by Feynman 
diagramsP^G2] Figure [3] shows the 2- and 3- loop diagrams contributing to 
In Z as an example. Each loop, which are considered to have no orientation, 
carries a factor of N. The contact interaction is represented by a dashed 
line for convenience, each carrying a weight oc g. In the equivalent space- 
time approach, Feynman diagrams are understood as spacetime diagrams. 
Each loop of a Feynman diagram is expressed in a hopping expansion as 
in the noninteracting theory, with the interaction standing for an inter- 
section between different worldlines or between different parts of a single 
worldline (self-intersection). In addition to the ones indicated, a worldline 
contributing to a given Feynman diagram may have further intersections, 
but these are treated as in the noninteracting case, i.e., as nonexisting. Put 
differently, apart from the intersections explicitly indicated, the particle 
trajectories still behave as phantom paths. 

Picturing the topology of closed worldlines, the first Feynman diagram 
(oc N) of the 2-loop contributions (oc g 2 ) in Fig. [3] represents single closed 
particle trajectories of arbitrary length and shape with one self-intersection 
carrying a factor of g. The second Feynman diagram (oc N 2 ) represents two 
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Fig. 3. Two- and 3-loop Feynman diagrams contributing to InZ of the \ip\ 4 theory. 

intersecting closed worldlincs. The 3-loop contributions (oc g 3 ) represent 
single closed worldlines with two self-intersections (1st and 2nd diagram 
oc N), two closed worldlines intersecting twice (3rd diagram oc N 2 ), two 
closed worldlines intersecting once with one of the two having in addition 
one self-intersection (4th diagram oc N 2 ), and finally three closed worldlincs 
(last diagram oc iV 3 ) connected by two intersections. 

Similarly, the nth order contributions to InZ, corresponding to the (n — 
l)th term (oc g n ~ 1 ) in the Taylor expansion (|6"6")) . all involve a total of 
ri—l intersections (including self-intersections) and up to n separate loops. 
The resulting loop gas is similar to that of the free theory, but the random 
loops are no longer phantom loops as the intersections explicitly indicated 
carry a weight oc g. In the path integral formulation, this is reflected by an 
additional steric repulsion in the worldline action, 

^ (i) =E [|^(4) + m 2 ] 

/ ds' r ds' r , S [x r (s' r ) - x r '(s' r ,)} , (67) 
__, Jo Jo 



9 
"6 



r,r' • 



where the sums J2 r anc ^ Sr' extend over the / particle trajectories assumed 
to be present. The net result is that (self-)intersections are suppressed and 
the random loops tend to become mutually and self-avoiding, the more so, 
the larger the coupling constant g becomes. Turning on the interaction thus 
changes the typical form of the random loops and in particular their fractal 
structure. 



To derive Eq. (|67|) . we follow SymanzilPfil and write the interaction 
term in the action ([2]) in the continuum limit as a functional integral over 
an auxiliary field a 



eXp (~4! / ddx(fi4 ^) = / D(7CX P - J d d x (^a 2 - iatp 2 



(68) 
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up to an irrelevant prefactor. The partition function JT]) then becomes in 
the continuum limit 



d d x 



1/^ -.2 m 2 

■^PnV) +—¥> -ur<P~ 



Dip Der exp 

J Dcr Vet- N/2 (-fig + m 2 - 21a) exp [~ J ^ 



oo 1 / N \l I 



1=0 



n 



n 



Day(0 



exp 



.9 



dW 



ex p( _ E/ r d4{ii 2 (4)+m 2 -2ia[ a ; r ( S ;)]} 



(69) 



where the last equality follows from the previous result (|64p . The auxiliary 
field (j is a function of the coordinates along the particle trajectories, a = 
a [x r (s' r )). A simple Gaussian integration yields 



Dcr exp | -- J d d xa 2 + 2\^ J ds' r a[x r (s' r )] 



exp 



- y 

ft Z-^i 



ds' r / ds' r ,6[x(s' r ) -£«,)] ^ (70) 
o Jo 



again up to an irrelevant prefactor. Using this in the last expression for Z, 
we obtain the spacetime representation of the partition function of the ip 4 
theory 



E 

1=0 



1 f N 
l\ 



i i 

n 



ds T 



Dx r (sJ,) 



-w 



(0 



(71) 



with the classical action of the point particles tracing out their world- 
lines given in Eq. (|67|) . The Ith term in this sum describes a tangle of I 
loops. Since the number of loops can be arbitrary large, Eq. (f7Tj) denotes 
the partition function of a grand canonical ensemble of fluctuation loops, 
each of arbitrary length and shape, with steric repulsion. 

1071 

Surprisingly, as was first pointed out by Balian and Toulousep- iJ the 
noninteracting theory can be recovered from the interacting theory by tak- 
ing the number of field components N — > —2. Consider a contribution to 
the partition function involving the first vertex in Fig. 01 Assume that leg 
1 is connected through the rest of the Feynman diagram to leg 2. Both 
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legs carry the same spin index a, which after taking the sum yields a fac- 
tor of N. Legs 3 and 4 are necessarily also connected through the rest of 
the Feynman diagram. They carry the same spin index /3, say, which after 
taking the sum also yields a factor of N. Replacing the first vertex with 
the second or third vertex in Fig. |4] yields identical contributions save for a 
change in topology regarding the routings of the spin indices. Specifically, 
the graph with the first vertex involves two different sets of spin indices: 
a and /3, whereas the graphs with the second or third vertex involve only 
a single set of spin indices as a — (3. Since the first graph carries an extra 
factor of N = — 2 (the loop fugacity) after the sum over the spin index is 
taken, the three graphs cancel. All other graphs involving vertices cancel 
in the same fashion three by three. What remains are phantom loops of a 
noninteracting theory. 





Fig. 4. The representation of the vertex of the tp theory by a dashed line makes the 
routings of the field components, labeled by a and /3, explicit. 



2.4. O(N) Spin Model 

Whereas the free theory discussed in Sec. 12.11 corresponds to taking the 
coupling constant A — > in the action ©, the O(N) spin model corresponds 
to taking the opposite limit, A — > oo. As for A —* 0, the theory becomes 
quadratic also in this limit of infinite repulsion. The theory is nevertheless 
nontrivial because of the fixed-length constraint 0^ = N at each lattice site. 
Instead of considering the conventional Boltzmann weight exp(K<f> n ■ </>«<), 
often a simplified representative of the 0(N) universality class is studied, 
obtained by truncating that weightPHl 

Z = Ti H {l + Kct> n -(j> n .). (72) 

(n,n ; ) 

The product is restricted to nearest neighbor pairs. The remaining factor in 
the on-site measure (fT0|) becomes trivial in this limit and has been ignored. 
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The main difference with the original spin model is that in the truncated 
model bonds cannot be multiply occupied. It is generally accepted that this 
simplification does not change the universal properties of the theory. 

The hopping expansion is usually referred to the high-temperature ex- 
pansion in the context of spin models. Also for the truncated model, it cor- 
responds to expanding the partition function in terms of the parameter K. 
The contributions to Z can again be visualized by graphs on the latticePSl 
As will be detailed below for the Ising model, only closed graphs with an 
even number of occupied bonds connecting each vertex yield nonzero contri- 
butions. The virtue of the truncated model is that bonds cannot be multiply 
occupied. In contrast, no restriction on multiple occupancies exists in the 
representation with the usual Boltzmann weight. When formulated on a 
two-dimensional honeycomb lattice, which has coordination number = 3, 
the high-temperature graphs of the truncated model are automatically mu- 
tually and self- avoiding, a given lattice site is either empty or has one 
bond entering and one leaving the site (see Fig. 




Fig. 5. Closed high-temperature graphs on a honeycomb lattice automatically form 
mutually and self-avoiding loops. 

The partition function can then be written geometrically as a sum over 
all possible mutually and self-avoiding (MSA) loopsp^l 

z=J2 NlRb > ( 73 ) 

MSA 
loops 

with b and I the number of occupied bonds and separate loops forming 
the graph. Equation (|73[) constitutes the high-temperature representation 
of the partition function. As for the noninteracting field theory ([63]) . each 
occupied bond carries a weight K (bond fugacity) and each loop carries a 
weight N (loop fugacity). Whereas the random loops of the noninteracting 
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field theory are phantom loops that can freely intersect and share bonds, 
here they are mutually and self- avoiding. As before, they physically denote 
the worldlines of the particles described by the quantum field theory. 
The O(iV) spin-spin correlation function 




is represented diagrammatically by a modified partition function, obtained 
by requiring that the two sites n and n' are connected by an open high- 

H <>| 

temperature graph. 1 ^ On a honeycomb lattice, the scaling part of the cor- 
relation function is given by the connected graphs 

~ E K " = ls a , p y Ez b (n,n')K b , (75) 

graphs b 

where Zb(n,n') is the number of (open) mutually and self-avoiding graphs 
along b bonds connecting the lattice sites n and n'. Strictly speaking, 
(finfin' ) < \&a,P Sgraphs ^ b as tne cancellation of the disconnected graphs 
in the numerator and Z in the denominator, required for an equality, is not 
complete: For a given open graph, certain loops present in Z are forbidden 
in the modified partition function as they would intersect the open graph, 
or occupy bonds belonging to it, which is not allowed. In other words, the 
presence of an open graph influences the loop gas and vice versa. Since 
each loop carries a factor N, the loop gas is absent in the limit N — > 0, 
pioneered by de Gennesff21 The inequality in Eq. |75|) becomes an equality 
in this limit, and the open graphs become ordinary self-avoiding random 
walks on a honeycomb lattice. These walks physically describe polymers in 
good solvents at sufficiently high temperatures so that the van der Waals 
attraction between monomers is irrelevant. For the noninteracting theory, 
where the worldlines are phantom trajectories, the loops cancel, so that the 
corresponding correlation function is exactly given by the open trajectories 
connecting the endpoints [see Eq. (fT6|) ] . 

For a d-dimensional hypercubic lattice, the coordination number is 
z = 2d, so that intersections are now possible. Although the bonds of the 
truncated model (|72|1 can by construction still not be multiply occupied, 
the high-temperature graphs are no longer mutually and self-avoiding. Af- 
ter rescaling the bond fugacity K — > K/N, the partition function (|72| can 
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be written aJ^Sl 




) 



ft 



Z = V T V T-" V : 



(76) 



loops 



where m 2 , • • • , rri2d indicate the number of intersections of 2k bonds (k = 
1, • • • , d) forming a given t ang le of b bonds and I loops. An intersection of 
2k bonds carries the weiglrF^ 



Due to the rescaling of the bond fugacity, = 1 for all N. In the limit 
N — > 0, only the weight V2 survives as V2k — > for k > 1, implying that 
even on a hypercubic lattice, intersections are completely suppressed and 
the high-temperature graphs reduce to self-avoiding walks. 

In the opposite limit, N — > oo, all the vertex weights become unity 
V2k - > lj and the partition function (|76|) takes the noninteracting form (|63p 
after undoing the rescaling of the bond fugacity. The only difference with 
the noninteracting theory is that for a given number of set bonds 6, only the 
graph with the highest number of loops needs to be included as it dominates 
all other graphs. This limit corresponds to the so-called spherical modelP^ 

The loop representation (|76p. featuring colorless loops, obtains by re- 
solving a vertex where 2k bonds of an high-temperature graph meet into 
all possible routings of the paths (see Figs. [5] and [7]) • The number of rout- 



Fig. 6. A vertex with four bonds (k = 2) resolved into three possible routings of the 
paths. 

ings equals the number of ways 2k objects can be divided into k distinct 
pairs, i.e., (2k — 1)!!. Each of the routings carries the same weight. A given 
high-temperature graph is thus resolved into loop configurations, and the 
partition function ((76|) describes a loop gas, as does Eq. (|73l) . The loops 
physically represent worldlines that can now intersect. Although loops on a 
honeycomb lattice cannot do so, for both the honeycomb and square lattice 
the same large-scale behavior is obtained. This implies that the operator 



N k T(N/2) N k 



(77) 



V2k = 



2 k T(k + N/2) N(N + 2) ■■■(AT + 2fc - 2) ' 




+ 
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Fig. 7. A graph (top panel) involving two vertices with four bonds (k = 2) is resolved 
into 9 different loop configurations according to the rule displayed in Fig. \E\ 



which introduces intersections is an irrelevant perturbation at the critical 
point SH 

It is worth pointing out that in the high-temperature representation, 
the observation that N — —2 describes phantom loops is far from obvious. 
On a square lattice, for example, a graph of b bonds and one vertex where 
four bonds meet contributes the following term to the partition function, 



lim Vi N(N + 2) I ^ 

\N 



(-2) 2 ( -? 



(78) 



which is finite because the weight factor given in Eq. (|77|) diverges in the 
limit N — > — 2. The diagrams no longer cancel three by three as they did 
in Sec. O 
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2.5. Ising Model 

As an explicit example, we consider the simplest 0(A) model, viz. the Ising 
model (N — 1) on a square lattice. The model is defined by the Hamiltonian 

H = -jJ2<t>n(l>n>, (79) 
(ran') 

where the spins can take only two values <p n — ±1, and the interaction 
parameter J is taken to be unity. The standard partition function reads 

Z = Tr c- fJH = Tr JJ e"*»*»' , (80) 

(ran') 

with the trace Tr denoting the sum over all possible spin configurations. 
In the context of statistical physics, (3 denotes the inverse temperature 
[3 = l/fceT, with fee Boltzmann's constant. In the present context, we 
regard (3 as coupling constant. Because of the identity 

e PtnK> = CO sh/? (1 + Kfafa) , (81) 

with K = tanh/3, the standard representation can be rewritten in the trun- 
cated form (172[) as: 

Z = (cosh/3) 2Ar Tr JJ (1 + (82) 

{nn') 

up to an irrelevant prefactor. Here, 27V denotes the total number of bonds 
on a square lattice with N sites and periodic boundary conditions. The right 
hand contains a total of 2N factors, one for each bond, and therefore a total 
of 2 2N terms when multiplied out. Each term can be visualized by a graph 
along the bonds on the lattice, with each bond representing a pair of adja- 
cent spins. Because of the trace in the definition of the partition function, 
a graph having a loose end does not contribute to Z. Indeed, let the spin 
at the loose end be 4> n then the sum over its possible configurations yields 
zero, X^„=±i 0« = 0. As a result, only closed graphs contribute. Similarly, 
only those closed graphs with an even number of bonds, i.e., two or four for 
a square lattice, connecting each vertex of the graph yield nonzero contri- 
butions. Such a vertex contributes a factor of two as J2<p =±i 4*™ = 2 for 
m even. The high-temperature representation of the Ising model therefore 
involves a sum over all possible closed graphs that can be drawn on the 
lattice: 

Z=(coshl3) 2N 2 N Y, K "- ( 83 ) 

closed 
graphs 
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A graph of b bonds can be built from various disconnected closed sub- 
graphs. The high-temperature representation is to be compared with the 
loop representation (|76p . which reduces to 

loops 

for N = 1 on a square lattice, with m denoting the number of vertices 
with four bonds. Remembering that the loop representation is obtained by 
resolving each of these vertices according to the rule given in Fig. [5J we 
see that apart from an irrelevant prefactor the two expressions coincide. 
The Ising model (N = 1) is special as the loop fugacity is unity, so that 
the contribution of a closed graph depends only on the number of bonds it 
occupies. 

The correlation function of the Ising model, 

(mm') 

= |(cosh/3) 2A, Tr0„^ Yl (\ + K^ m ,). (85) 

(mm') 

can similarly be visualized by graphs on the lattice. Each graph now con- 
tains a subgraph connecting the lattice sites n and n' . All other subgraphs 
must be closed and, except for the endpoints, all vertices must be even as 
before. In terms of these graphs, 

{KM = | (cosh Z?) 2 ^ Y, K "- ( 8fi ) 

graphs 

2.6. Summary 

In the spacetime approach to quantum field theory on a lattice, which is 
geometric in nature, the partition function and the correlation function are 
calculated in a hopping expansion, thus reducing the problem to counting 
paths on the lattice. Physically, the paths represent the worldlines traced 
out by the particles when hopping from one lattice site, representing a 
spacetime cell, to an adjacent site. The noninteracting lattice field the- 
ory features phantom worldlines that can freely intersect and share bonds. 
When the self-interaction is turned on, by taking the coupling constant 
A > 0, intersections carry an energy penalty. This has the effect to suppress 
such configurations and as a result the particle trajectories tend to avoid 
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themselves. In the limit of infinite repulsion A — > oo, the (ft 4 lattice field 
theory reduces to the O(N) spin model. The hopping expansion, known in 
this context as the high-temperature expansion, takes a particular simple 
form when the truncated O(N) model (|72p is implemented on a honeycomb 
lattice. The worldlines are then by construction mutually and self-avoiding. 
Since the large-scale behavior of this representative of the O(N) universal- 
ity class on the honeycomb lattice is the same as for the generic <p A field 
theory defined on a square lattice, both multiply occupancy of bonds and 
intersections are irrelevant. At large scales, both loop gases have the same 
fractal structure. Even when considering weak repulsion by taking A small, 
after coarse-graining the fractal structure of the loops appears the same as 
in the limit A — ► oo. This limit therefore governs the large-scale behavior. 

3. Critical Properties 

In this section, the critical properties of the 0(JV) universality class are dis- 
cussed from the perspective of worldlines. In particular, the fractal structure 
of the loop gas featuring in the spacetime approach to fluctuating fields close 
to the critical point is studied. It is shown how the critical exponents char- 
acterizing the phase transition can be obtained from the fractal structure of 
these geometrical objects. The discussion is geared after percolation theory, 
which is purely geometric in natureP^ 

3.1. Fractal Structure 

Consider the O(N) theory close to the point where it undergoes a contin- 
uous phase transition to the ferromagnetic state, where the spins prefer to 
point in one direction in spin space. This phase transition is associated with 
the spontaneous breaking of the O(N) symmetry to the subgroup 0(N — 1) 
which includes rotations about the preferred spin direction. The logarithm 
of the partition function of the interacting theory can in the vicinity of this 
critical point be written in a form similar to the one found in Eq. (|54|) for 
the noninteracting theory as 



b 

The worldline loop distribution £{, asymptotically splits into two factors as 
it did for the phantom loops in Eq. (|55j) : 




(87) 



b 



(88) 
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with t an exponent characterizing the interacting theory. The value r = 
d/2 + 1 found for phantom loops is typical for a noninteracting theory where 
the loops are simple closed Brownian trajectories. Formally, Eq. (|87|) is the 
same as for clusters in percolation theory close to the percolation threshold. 
As in that context J^r can be related to the fractal dimension of the loops. 
Consider the square radius of gyration B% of a loop of b steps 

R i = \ E>* - n ? = i E fa ~ n ^ 2 ' ( 89 ) 

i—l i,j=l 

with rii denoting the lattice sites visited by the particle while tracing out 
its worldline on the spacetime lattice, and n = (1/6) 53i=i n * the cen ter of 
mass of the loop. Asymptotically, the average (R%) scales with the number 
of steps b taken as 

(Rl)~b 2 / D , (90) 

which defines the Hausdorff, or fractal dimension D. This average yields the 
typical linear size of a loop of b steps. Specifically, such a loop is distributed 
over a volume of typical linear size (i?^) 1 ^ 2 . Given the definition of the loop 
distribution as the average number (per lattice site) of loops of b steps 
present, it follows that bl^ is the probability that a randomly chosen lattice 
site belongs to such a loop, and 

b£ b ~ 1/ (R 2 b ) d/2 , (91) 
with d the dimension of the spacetime box. This leads to the relation 

r = - + 1, (92) 

which is formally the same as in percolation theory!^ 

3.2. Loop Proliferation 

The Boltzmann factor in the loop distribution (|88[) exponentially suppresses 
long loops as long as the line tension 8 is finite. Upon approaching the 
critical point, the line tension vanishes as a power law, 

~ (K c - Kf'\ (93) 

with a a second exponent characterizing the loop distribution. When this 
happen s, th e loops can grow arbitrarily long at no cost and the entire system 
becomes^ "pierced through and through with" worldlines. It is through this 
loop pro/i/erataonloop!proliferatiorl^lthat a phase transition manifests itself 
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in the spacetime approach. The phenomenon is completely analogous to the 
sudden appearance at the percolation threshold of a cluster spanning the 
infinite lattice in percolation phenomena. 

In the field theoretic approach, the phase transition manifests itself 
through a nonzero vacuum expectation value of the field 0", where the 
component a indicates the preferred spin direction spontaneously chosen 
by the system. Since all possible directions are equivalent, the system can 
choose any of these. On a finite lattice, the system can easily change its 
preferred direction. This is in particular the case for Monte Carlo simula- 
tions using a nonlocal update, such as the Swendsen- Wang 331 or WbmSSl 
algorithm, in which not individual spins are updated, but entire clusters. 
For this re ason , the following order parameter is frequently used in lattice 
simulations!^ 




With this choice of the order parameter, all possible directions are treated 
equally. In the normal phase, the expectation value of this operator vanishes. 
At the critical point, a condensate forms and a nonzero expectation value 
develops, 

(O) ± 0, (95) 

typical for the ferromagnetic state with the spins pointing in a preferred 
direction in spin space. Those particles tracing out long worldlines all belong 
to the condensate. 

The loop distribution ([55]) is related to the number of paths Zbifl) of b 
steps returning to a site adjacent to the initial position through 

4 = \z b {a)K b . (96) 
o 

Because z b (a) refers to closed worldlines, the factor 1/6 is included to pre- 
vent overcounting as a given loop can be traced out starting at any lattice 
site along that loop. It should be noted that in Eq. (|96|) rooted loops are 
considered, i.e., all loops run through an arbitrary but fixed lattice site n. 
The relation (l96|) implies that Zb(a) scales as 

z b {a)K b ~ b- d/D e~ eb . (97) 

Similarly, the number of paths z b of b steps starting at n and ending at an 
arbitrary lattice site defined in Eq. (|40|) scales as 

z b K b ~ b^ D e' eb 1 (98) 
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with i? a new exponent, characterizing the number of ways an open path 
of b steps can be embedded in the lattice. It is worth pointing out that in 
contrast to the loop exponents r and a, this third exponent refers to open 
paths. Its value is expected to be positive because the number of possible 
rooted open paths with no constraint on their endpoints increases with the 
number of steps b. This is in contrast to the loop distribution (|88[) . where 
the algebraic factor decreases with increasing b, reflecting the increasing 
difficulty for a path to close. According to Eq. (f3"9")> . the ratio of Zb(n,n') 
and Zb defines the probability Pb(n,n') for a particle to move from the 
lattice site n on the spacetime lattice to n' in b steps. On general grounds, 
it scales at criticality as: 

P b (n, n') ~ b- d /° P (\n - , (99) 

with P a scaling function. This generalizes the Gaussian distribution (|44|) 
found for the noninteracting theory. Consistency of the three scaling for- 
mulas (l97|) . (|98|) . and (|99p. requires, as was first shown by McKenzie and 



MooreS-^for self-avoiding random walks, that the scaling function P(t) must 
vanish for t — ► as a power law, 

P(t) ~ t* (100) 

with the same exponent governing the asymptotic behavior (|98p of the num- 
ber Zb of open paths at the critical point (9 = 0). 

3.3. Critical Exponents 

As in percolation theory, the loop exponent a is related to the exponent v, 
which specifies how the correlation length £ diverges at the critical point, 

t~\K-K c \- v . (101) 

Physically, £ indicates the typical length scale in the system. The radius of 
gyration-a second typical length scale-can be written in terms of £ as 

(R b )=tR(b9), (102) 

where R is a scaling function and 9 is the same parameter as in the loop 
distribution . From the asymptotic behavior (|9"0|) , the divergence of the 
correlation length, and the vanishing of the line tension 9 as K c is 
approached, the relation^] 

, = -L, (103) 
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or 

v = T —± (104) 
da 

follows. These relations are formally identical to those in percolation theory. 

The rest of the thermal critical exponents can be equally well expressed 
in terms of the loop exponents and $. To derive these expressions, we write 
the correlation function in terms of zi,(n,n') as on the honeycomb lattice 
[see Eq. ([75]) ] 



'urn ^ 

b 



When evaluated at criticality, where 

/^J-WI"-"'!^'. (we) 

Eq. (fT05j) gives 

r) = 2-D-$. (107) 

This relation, with <d defined by Eq. ()L00|1 . was recently proposed for the 
three-dimensional XY, i.e., 0(2) model in Ref. [48]. This model describes 
the superfluid phase transition in 4 He. 

Finally, using the definition of the magnetic susceptibility %, 



X 



which diverges as x ~ I K ~ Kc | 7 , we find 

7=~(l+£). d09) 

The explicit expressions given for is, ry, and 7 satisfy Fisher's scaling relation, 
j/u = 2 — rj, and 

- = D + d. (110) 
v 

With v = (t — I)/ 'da, Eq. {53]) yields the scaling relation dv = 2 — a, 
where a determines the scaling behavior of the free energy close to the 
critical point, 

lnZ/fi- \K -K c \ 2 ~ a . (Ill) 

The expressions for the other exponents follow by using the remaining scal- 
ing relations. We thus have shown that all the thermal critical exponents are 
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determined by the configurational entropy exponents for closed and open 
particle trajectories, r and <&, respectively, and by a. 

It is worth noting that the thermal critical exponents depend on only 
two independent variables, viz. D + i? and aD. Their significance in field 
theory is that they determine the anomalous scaling dimensions of the <p 
and if 2 fields 

d v = \{d-D-<&), d (p 2=d-aD, (112) 

respectively. 

3.4. Self-Avoiding Random Walks 

Before applying the results of the previous section to the two-dimensional 
O(N) spin model with arbitrary — 2 < N < 2, let us first verify that these 
general results are consistent with the results known in the polymer limit 
N -> fjfSl where, as we saw below Eq. ([77]). the high-temperature graphs 
reduce to self-avoiding walks. Through the exact enumeration and analysis 
of the number of self-avoiding loops on a square lattice up to length 110, 
the value r = 5/2 for the loop exponent r has been established to very high 
precision!^ According to Eq. ([52]) . this corresponds to the fractal dimension 
D = 4/3-a value that has been independently established to high precision 
in that same study by determining the average square radius of gyration 
{R 2 } of the loops. In a related studyf^where the number Zb of self-avoiding 
walks on a square lattice up to length 71 has been enumerated and analyzed, 
the value of the exponent characterizing the open trajectories has been 
established. Specifically, the value d/D — 11/32 was found to high precision, 
yielding = 11/24. 

In most studies on self-avoiding walks, the fractal dimension of the 
walks is simply equated to the inverse correlation length. As follows from 
Eq. (|103p . this corresponds to setting a = 1. A closer inspection of dc 
Gennes' original work on the N — > limit reveals that there indeed a = 1. 

To sum up, the two configurational entropy exponents r (loops) and & 
(open paths), and a for a self- avoiding walk in two dimensions are given by 

r = 5/2, i?= 11/24, a=\. (113) 

With these values, all the critical exponents of the two-dimensional 0(N — > 
0) model can be found through the relations of the previous subsection and 
standard scaling relations. 
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3.5. O(N) Models 

The class of two-dimensional O(N) models for —2 < N < 2 can be 
parametrized as 



leading to r = (6 + R)/(2 + R). The observable whose scaling dimension is 
given by this fractal dimension consists of two spins in different spin states 
placed at the same site: 0"</>^, with a ^ /^PSl It physically measures the 
tendency of spins to align, as becomes intuitively clear by remembering that 
the spins connected by a high-temperature graph are all in the same spin 
state. 

While the fractal dimension of the high-temperature loops is known 
exactly in two dimensions, the value of the entropy exponent for open 
paths has to our knowledge not been directly established for N ^ 0. In 
two dimensions, it can be established indirectly^! by considering the Potts 
model-a closely related spin model. 

That model was reformulated by Fortuin and KasteleyrJ^ in a purely 
geometric fashion in terms of suitably defined spin clusters. The thermal 
phase transition is marked by a proliferation of these clusters, which more- 
over have the thermal critical exponents encoded in their fractal structure. 
Percolation observables such as the average cluster size and the percolation 
strength, defined as the probability that a randomly chosen site belongs 
to the percolating cluster, serve as improved estimators for the magnetic 
susceptibility and the magnetization, respectively. In other words, the ther- 
mal critical exponents of the Potts model can be determined by studying 
the fractal structure of these so-called Fortuin-Kasteleyn clusters, with the 
exponents defined as in percolation theory. In two dimensions, the fractal 
dimension Dq of the Fortuin-Kasteleyn clusters is known exactlj®l 



where the Q-state Potts models are parametrized in terms of the same 
parameter R parametrizing the 0(JV) models as 




(114) 






(116) 




(117) 
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with | < R < 1, so that < Q < 4. The Q -> Potts model (R = |) 
describes standard percolation theory. Surprisingly, the external perimeters 
of the Fortuin-Kasteleyn clusters have the same fractal dimension D as the 
O(N) loopsP^For example, the external perimeters of clusters in standard 
percolation theory have the dimension | as have self-avoiding random walks 
(k = |). It can be shown that the entropy exponent $ for open trajectories 
and the two Potts dimensions satisfy the sum rule^l 

d = 2D C - D - i9, (118) 

with d = 2 the dimension of the spacetime box. With D and Dq known 
exactly, the expression for {) follows. 

Finally, the expression for the exponent a can be inferred from the values 
a = 1 for N -> (re = |) and a = for JV = 2 (k = 1). This last value 
reflects the special status of the XY model, undergoing a phase transition 
of the Berez inskii- Kosterlitz-Thouless type characterized by algebraic long- 
range order.ESES] Assuming, as is the case for Fortuin-Kasteleyn clusters, 
that both r and a have a similar dependence on k, we deduce that a — 
8(1-k)/(2 + S). 

To sum up this subsection, the two configurational entropy exponents 
r (loops) and § (open worldlines), and a for the O(N) model with — 2 < 
N < 2 are given in terms of R by 

6 + k„ k31 1 — R , 

t = ——, $ = -! + - + — , cr = 8 . (119) 

2 + R 7 2 4k' 2 + R V ' 

Table [T] collects the resulting values for the various two-dimensional O(N) 
models together with the corresponding thermal critical exponents. It is 
worth pointing out that for N — — 2 indeed Gaussian exponents are ob- 
tained, but that the fractal dimension of the high-temperature graphs is 
not given by that of a Brownian random walk, which has D = 2. Also note 
that with increasing N, the exponent d decreases, while at the same time 
the fractal dimension D increases. The ratio fi/D appearing in the scaling 
relation (|98|) thus decreases with increasing N, implying that the number of 
possible rooted open paths with no constraint on their endpoints increases 
less rapidly with the number of steps b for larger N. 

3.6. Summary 

It was shown that the worldlines or high-temperature graphs of the O(N) 
theory proliferate right at the thermal critical point, and that the entire set 
of critical exponents can be retrieved from the fractal structure of these line 
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Table 1. Critical exponents together with the exponents characterizing the 
high-tcmperature graphs (worldlines) of various two-dimensional O(N) models. 



Model 


N 


K 


a 


P 


7 




V 


D 


T 


■3 


(7 


Gaussian 


-2 


1 

2 


1 





1 





1 

2 


5 
4 


13 
5 


3 
4 


8 
5 


SAW 





2 


i 


5 


43 


5 


3 


4 


e 


11 


1 


3 


2 


04 


:V2 


24 


4 


3 


2 


24 


Ising 


1 


3 

4 





1 
8 


7 

4 


1 

4 


1 


11 

8 


27 
11 


3 
8 


8 
11 


XY 


2 


1 


— oo 


OO 


CO 


1 

4 


OO 


3 
2 


7 
3 


1 

4 






objects. In this way, a purely geometric description of the phase transition 
in these systems was arrived at. 

4. Monte Carlo Simulations 

We next discuss a Monte Carlo study of the Ising model on a square lattice 
we carried oulPQI to support the above findings (for a summary of that 
study, see Ref. [57] )• We consider not the original spin formulation of the 
model, but the high-temperature representation detailed in Sec. l2.5l instead. 
and examine whether the fractal structure of the high-temperature graphs 
indeed encodes the thermal critical behavior. We focus exclusively on the 
graph configurations and use percolation observables to study the fractal 
properties of the graphs. 

Traditionally, high-temperature expansions are carried out exactly up to 
a given order by enumerating all possible ways graphs up to that order can 
be embedded in the lattice. 1 ^ We take a different approach and study the 
high-temperature graphs by means of Monte Carlo methods. The Metropo- 
lis algorithrrP^l W e use involves a local update and for that reason suffers 
from critical slowing down close the the transition temperature. We have 
intensionally chosen the local update as it allows for a particular clean im- 
plementation of our ideas. The price we pay for this is that we cannot study 
too large spacetime boxes. The largest lattice we can reasonably study is of 
linear size L — 256. 

4.1. Plaquette Update 

The update algorithm we use in our purely geometric approach of the Ising 
model directly generates the high-temperature graphs that contribute to the 
partition function. Specifically, the closed graphs are generated by means 
of a Metropolis plaquette algorithm. A proposed plaquette update resulting 
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in V 



occupied bonds is accepted with probability 



Pht = 



{ 



K' 
1 



if b 1 > b 
else 



(120) 



where b denotes the number of occupied bonds before the update!^ In the 
entire temperature range < (3 < oo, K — tanh/3 < 1. Reflecting the Z2 
symmetry of the model, all bonds of an accepted plaquette are changed, i.e., 
those that were occupied become unoccupied and vice versa (see Fig. [8]) . 
With &□ denoting the number of bonds on the proposed plaquette already 



Fig. 8. Plaquette update: An existing high-temperature graph with the plaquette pro- 
posed for updating indicated by the broken square (left panel) and the new graph after 
the proposal is accepted (right panel). 

occupied, it follows that 



The prescription (|120p immediately follows from the high-temperature 
representation (|83[) of the Ising model. In equilibrium, the probability distri- 
bution P(G) for a given graph configuration G involving b occupied bonds 
reads 



Such a configuration can be reached after t + 1 iterations from the config- 
uration present after t iterations in the following fashion 



P(G, t + l) = P(G, t) + J2 [ p ( G '> t)W(G' -»■ G) - P(G, t)W(G -> G')] , 



where W(G — ► G 1 ) is the probability for the system to move from the graph 
configuration G with b occupied bonds to the graph configuration G' with 
b' occupied bonds. In equilibrium, P(G,t + 1) = P(G,t) = P(G), and the 
system satisfies detailed balance 



b' = b + 4-2b a . 



(121) 



P(G) = -{cosh(3) 2N 2 N K b . 



(122) 



G' 



(123) 



P(G')W(G' -»G) = P(G)W(G -* G'), 



(124) 



February 6, 2008 22:45 



WSPC/Trim Size: 9in x 6in for Review Volume 



3ding 



36 W. Janke and A. M. J. Schakel 

or 

W(G -> gO _ P(G") _ K b ' 

W(G' -> G) ~ P[G) ~ K b ' (1 5) 

To maximize the acceptance rate of proposed updates, the largest of the two 
transition probabilities W[G — > G') and W(G' — > G) appearing in the ratio 
should be given the largest possible value, i.e., one. Thus, if the number of 
bonds b' in the proposed configuration is larger than the number of bonds 
b in the present configuration, so that K b /K b < 1, then W(G' — > G) = 1 
and W(G — > G 1 ) = pwr- On the other hand, if b' < b, the proposed config- 
uration carries a larger weight than the present one, and will be accepted 
unconditionally. 

By taking plaquettes, i.e., elementary loops on the spacetime lattice as 
building blocks, the resulting high-temperature graphs are automatically 
closed (see Fig. |9] for snapshots at three different temperatures) . 

4.2. Numerical Results 

The data was collected on lattices varying in linear size from L — 16 to 256 
in 3.3 x 10 5 Monte Carlo sweeps of the lattice close to the critical point and 
1.1 x 10 5 sweeps outside the critical region. About 10% of the sweeps were 
used for equilibration. After each sweep, the resulting graph configuration 
was analyzed. Statistical errors were estimated by means of binning. 

We first determine whether the high-temperature graphs proliferate pre- 
cisely at the inverse critical temperature j3 c = ln(l + y/2) /2 = 0.440686 
Remember that the inverse temperature is related to the tuning parameter 
K via K = tanh j3. To this end, we measure the so-called spanning proba- 
bility Ps as a function of (3 for different lattice sizes. Giving the probability 
for the presence of a graph spanning the lattice, Ps tends to zero for small 
f3, while it tends to unity for large [3. This observable has no scaling dimen- 
sion and plays the role of the Binder cumulant in standard thermodynamic 
studies. When plotted as a function of (3, the curves Ps(/3) obtained for 
different lattice sizes should all intersect in one point. Being independent 
of the lattice size, this common point marks the proliferation temperature 
of the infinite system. Within the achieved accuracy, we observe that all 
the measured curves cross at the thermal critical point, implying that the 
high-temperature graphs loose their line tension and proliferate precisely 
at the Curie point (see Fig. [T0|. 

According to standard finite-size scaling, this observable obeys the seal- 
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Fig. 9. Snapshots of high-tcmperature graphs of the Ising model on a square lattice with 
linear size L = 128 and periodic boundary conditions in the normal phase at fi = 0.6/3 c 
(top left panel), at the critical point f3 = f3 c (top right panel), and in the condensed 
phase at j3 = 1.4/3 C (bottom panel). Because of its richer structure, the last snapshot is 
enlarged relative to the first two snapshots. 
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0.8 1 1.2 1.4 1.6 1.8 



0/0c 

Fig. 10. Probability Pg for the presence of a spanning graph as function of the inverse 
temperature /3 measured for lattice sizes L = 16, 32, 64, 128, 256. Within the achieved 
accuracy, the curves cross at the thermal critical point (3 = j3 c . 

ing 

P s (f3,L)=P s (L/0, (126) 

where Ps is a scaling lunction and £ the correlation length whose divergence 
at criticality is governed by the exponent v. The data gathered on lattices 
of different sizes is therefore a function not of (5 and L separately, but 
only of the ratio of the lattice size and the correlation length. All the data 
should therefore collapse onto a single curve when plotted as a function of 
(/3//? c - \)L X I V . With v given the Ising value v = 1, this is indeed what we 
find (see Fig. [TU). 

The fractal dimension of the high-temperature graphs is best determined 
by following standard percolation theory - and measuring the percolation 
strength , defined as the fraction of bonds in the largest graph, and the 
average graph size xg (see Figs. IT3l and fT2|). The percolation strength, 
which is finite in the low-temperature phase (/3 > (3 C ), vanishes as a power 
law when the percolation threshold is approached, 

P°~C9-/3c) /3a . (127) 

At the same time, the average cluster size diverges, 

Xg~ |/3-&r 7G - (128) 

Here, (3q and 7g are two percolation exponents which should not be con- 
fused with the thermal critical exponents. The ratios /3g/^ and 7g/^ are 
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Fig. 11. Raw data of Fig. UOI replotted as a function of (/3//3 c - l)L x /", with the Ising 
choice v = 1. The data collapse is satisfactory over the entire temperature range. 




Fig. 12. Percolation strength as a function of the inverse temperature j3 measured on 
lattices of linear size L = 16, 32, 64, 128, 256. Note the logarithmic scale on the vertical 
axis. 



expressed in terms of the entropy exponent r of the graph distribution as 
in percolation theorjH^l 



^ = d- ~ 2 



— = d- ~ T 



(129) 



v t — 1 v r — 1 ' 
with d the dimensionality of the lattice. According to TableQ] r = 27/11 = 
2.4546 . . . for the two-dimensional Ising model. Since D = d/(r — 1), the 
fractal dimension is related to these percolation exponents via 



D = d-^ = 1 -(d+^), 
v 2 V v J 



(130) 
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x 




1 1.5 2 2.5 3 



Fig. 13. Average cluster size xg> excluding the largest cluster, as a function of the 
inverse temperature measured on lattices of linear size L = 16, 32, 64, 128, 256. 

again as in percolation theory!^ Close to the percolation threshold, the 
percolation strength and average graph size obey the finite-size scaling laws 

P G = i -/W„pG (i/0j Xg = L 7g/^ Xg (L/£), (131) 

so that the exponent ratios and thus D can be determined by considering 
the system at criticality, where these scaling relations imply an algebraic 
dependence on the system size L. 




Fig. 14. Rescaled raw data of Fig. 1121 plotted as a function of (/5//3 c — VjL 1 /", with 
the Ising choice v = 1, and /3q = 0.626. The data collapse is satisfactory in the critical 
region as well as in the entire high-temperature phase j3 < (3 C - The flaring out of the 
data points for /3 > /3 C marks the end of the critical region. 
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We numerically arrived at the estimated^ p G = 0.626(7), 7g = 
0.740(4), leading to a = 0.732(6), r = 2.458(5) in perfect agreement with 
the exact values a = 8/11 = 0.7273 . . . , r = 27/11 = 2.4546 . . . , and the 

I r i I 

predicted fractal dimension^ D = 11/8 of the high-temperature graphs. 
The data were fitted over the range L = 16 — 128, using the least-squares 
Marquardt-Levenberg algorithm. With the obtained value for (3q, a satisfac- 
tory collapse of the percolation data is achieved (see Fig. IT4|) . The collapse 
of the average graph size data is less satisfactory. This is not untypical for 
this type of observable. A similar bad data collapse we observed for the 
average size of geometrical (as opposed to Fortuin-Kasteleyn) spin clusters 
in the two-dimensional Ising model. 

4.3. Summary 

The case study of the Ising model on a square lattice in the high- 
temperature representation showed that the fractal structure of the high- 
temperature graphs indeed encodes the thermal critical behavior. This rep- 
resentation therefore provides a purely geometric and equivalent description 
of the original spin formulation of the model. 

5. Further Applications 

The loop gas approach is by no means limited to the O(A) theory. The 
formalism is general and can be applied to a host of other phase transitions, 
not necessarily involving the proliferation of line objects. Phase transitions 
involving the proliferation of surfaces and hypersurfaces can be treaded 
similarly. In this section, a few further applications are briefly discussed. 

5.1. Higgs Model 

A straightforward extension of the O(A) theory follows by considering elec- 
trically charged particles. The continuum action describing the charged sys- 
tem reads 



S = / d d 



(132) 

where F^u = d^A v — d u A^ is the electromagnetic field strength, with A M 
the vector potential in d spacetime dimensions. The scalar field ip is now 
complex with A/2 complex, i.e., A real components, where A is even. The 
action is obtained from the neutral theory by minimal coupling the scalar 
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field to the gauge field A^, with e the electric charge. The Maxwell term 
provides the standard kinetic term for the gauge field, while the last term 
with parameter a fixes a Lorentz-invariant gauge. The theory (|132p is known 
as the Abelian Higgs model in high energy and as the Ginzburg-Landau 
theory in condensed matter physics. The partition function, whose general 
definition is given in Eq. JT]), of the charged theory involves in addition to 
the functional integral over the scalar field also one over the gauge field, 
i.e., 

Tr = J DtpDA^, (133) 

where the functional measure J D(p is the continuum limit of the discrete 
measure (j4|) defined on the spacetimc lattice, with an analogues definition 
for J DA M . 

To arrive at the spacetime description of this theory, we concentrate on 
the coupling to the gauge field by considering only the first two terms in 
the action, 

1 ... . . ,2 m 2 
-\{d IJ ,-ieA IJ )if\ 

This action describes otherwise free particles coupled to an electromagnetic 
background field, specified by the gauge field A^. The corresponding parti- 
tion function 

hxZ = lnDet _JV/2 - (<9 M - ieA^f + m 2 j (135) 

is readily written as a path integral by interpreting the operator 
[da — ieA^) + m 2 as the Hamiltonian operator of a nonrelativistic charged 
particle of mass M — | moving in a constant potential V — m? in d 
dimensions in the presence of an electromagnetic background field. The 
corresponding classical action is 

W = f da' {\± 2 {s') -iex^[x(s')]+m 2 } . (136) 
Jo 

According to the mnemonic discussed below Eq. (|60p . the spacetime rep- 
resentation of the partition function (|135|) is given in terms of this action 
as 

lnZ = yfi ^ j>Vx{s')e- Wo . (137) 

Since N is even here, the simplest Higgs model corresponds to the value 
N = 2, which in the context of condensed matter physics provides an 
effective description of ordinary superconductivity. 
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It follows from Eq. (|137[) that a particle trajectory T picks up an extra 
phase factor 

U(T) = e ie -/o d *' A P- < W x ( a ')] — gie / r ds„.A„(s)_ (138) 

This phase factor results in a Biot-Savart type of interaction between two 
line elements dx^ , dy v of charged loopsPl as can be illustrated by consid- 
ering the average of the phase factor (|138|) 



(U(T)) = J DA M U(r) exp j^d 
e 2 f 

~2 J ^ x ^^V vD ^ x ~v) 



4 » v 



1 

2^ 



exp 



where D^ v is the correlation function of the gauge field 

d d k 1 



(2ir) d k 2 



Ak-x 



(139) 



(140) 



as can be read off from the action (|132[) . This long range Biot-Savart inter- 
action between line elements of charged loops sets the Higgs model apart 
from the neutral O(N) theory, where the loops experience only a steric 
repulsion. 

The 1^1 4 term in the action (|132[) can be treated as in the O(N) theory, 
leading to a steric repulsion for charged loops. Pasting the pieces together, 
one arrives at the spa cetime representation of the partition function of the 
Higgs modeffSEEHEQ] 



BA U exp 



d d x 



i i 



9. 



— - (f Dx r (s' r 
s r J 



«e I n 

/=0 V / r=l 

with the Z-particle action 



(141) 



- y 



dsi 



ds' r , 5 [x r (s' r ) — X r '(s' r /)] 



(142) 



The partition function p41[) represents a grand canonical ensemble of fluc- 
tuating closed worldlines of arbitrary length and shape traced out by elec- 
trically charged point particles. Since the matter field ip is complex, the 
worldlines have an orientation now. 
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From this spacetime representation it follows that the logarithm of the 
partition function can again be cast in the general form (|87| , with the loop 
distribution now describing charged loops. Because of the additional long 
range Biot-Savart interaction, the fractal structure of these loops will be 
different from those featuring in the neutral theory, which experience just 
a steric repulsion. The physical picture remains unchanged, however. As in 
the neutral theory, only a few small loops are present in the normal phase 
due to the finite line tension 6. On entering the Higgs phase, characterized 
by a sign change in the mass term of the Higgs model, the line tension 
vanishes, and the charged loops proliferate as they can become arbitrary 
long at no cost. The vacuum becomes filled with "a spaghetti of tangled 



5.2. Bose-Einstein Condensation 

The spacetime approach was originally formulated to describe nonrelativis- 
tic quantum mechanics^ which fundamentally differs from its relativistic 
counterpart. A relativistic quantum particle roams the space dimensions as 
well as the time dimension, so that its worldlines are random walks in space- 
time. Such particle trajectories are parameterized by the Schwinger proper 
time parameter, i.e., by their arc length. A nonrelativistic particle, on the 
other hand, executes a random walk in space only. Its time coordinate is 
given by absolute Newtonian time, which elapses uniformly forward. In a 
nonrelativistic ensemble, all the particles execute a random walk in space 
synchronously in time, with all the worldlines parameterized by Newtonian 
time. That is, the Schwinger proper time parameter gets replaced by New- 
tonian time in the nonrelativistic limit. At the technical level, this can be 
seen by taking the nonrelativistic limit c — > oo in, for example, the correla- 
tion function (|37p of the relativistic theory. The integral over s can in this 
limit be approximated by the saddle point!^ 



showing that the arc length parameter s of relativistic worldlines indeed 
becomes replaced by Newtonian time r. This difference has profound im- 
plications. Because Newtonian time always elapses forward, nonrelativistic 
worldlines can no longer form simple loops on an infinite lattice. As a result, 
the number of particles present in the system is fixed. On a finite lattice, 
worldlines still have the possibility to close by wrapping around the lattice. 

Figure [15] shows an example of the motion in imaginary time of three 
particles in one space dimension. Up to now, the finite lattice sizes were 



loops"Pl 



s = r/2m, 



(143) 
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Fig. 15. Worldlines of three nonrclativistic particles on a square lattice, representing 
one space and the imaginary time dimension, moving synchronously from the bottom 
to the top of the lattice. The particles execute a Brownian random walk in x-direction. 
Periodic boundary conditions are chosen in the imaginary time direction. 



considered a mere approximation to systems of infinite extent, in space as 
well as in time. In this subsection, we instead assume that the imaginary 
time dimension is finite. With the time variable r restricted to the interval 
< r < H(3, and (3 the inverse temperature (3 — one describes the 

system at finite temperature. The zero-temperature limit is recovered by 
letting the imaginary time dimension become infinite. For bosons, periodic 
boundary conditions are to be imposed in the imaginary time dimension, 
implying that the configurations of an ensemble of nonrelativistic bosons 
at imaginary time and at h(3 are identical. (In contrast to the rest of the 
paper, we explicitly display factors of h in this subsection.) 

The partition function Zi of such an ensemble of I identical nonrel- 
ativistic bosons can be expressed in terms of the probability amplitude 
x(0), x(t)1 for a particle to move from x(0) to x(r) in imaginary time r 
Ml 



G 
as 



Zi = 



J2 [ n ddx rG[x 1 (0) ) x P(1) (S/3)]---G'[x i (0),xp (0 (?i/3)], 

P r=l 



(144) 



where the final configuration {xp( r )(?i/3)} at r — hp is a permutation P 
of the initial configuration {x r (0)} at r = 0. In contrast to the convention 
used up to now, here, d denotes the number of spatial dimensions. All 
possible permutations of the starting points are included as indicated by the 
sum J2p- I R Feynman's spacetime approach, the amplitude G[x(0),x(r)] 
is written as a sum over all possible particle trajectories connecting the two 
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spatial endpointJH 

G[x(0),x(fi/3)] = J Dx(r)exp (-- dr {yx 2 (r) + F[x(r)]}j , 

(145) 

where V[x(t)] is the potential experienced by the particle along its w orld line 
x(r). Expressed as a path integral, the partition function thus readJ^ 



r 1 / i 

Zt = e / n dx -( t ) ex p - * / drL(i) 

P •* r=l \ J° 



(146) 



with the Z-particle nonrelativistic Lagrangian 

£ (0 = £ + 5 J2 V[\Mr)-^(r)\] (147) 

and pair potential y(x r ,x r /) = T^(|x. r — x r /|). Figure H~51 gives an example of 
a configuration contributing to Z%. Since Xp^^fi/?) = xi(0), Xp( 2 )(fi/3) = 
X2(0), and Xp( 3 \(7i/3) = X3(0), the three particle trajectories form three 
separate loops around the time axis. This is typical for the high-temperature 
phase of a nonrelativistic ensemble of bosons, where most particles execute 
a random walk in space during the time interval < r < h(3 that returns 
to its initial position. It means that the particles, being distinguishable, 
behave classically. 

By mapping Feynman's spacetime approach onto the problem in classi- 



cal statistical mechanics of chains of beads connected by springs^ to nu- 
merically evaluate the path integral with the help of the Metropolis Monte 
Carlo method, Ceperley and Pollockj^SI and others showed that a powerful 
computational tool emerges that can even accurately describe a strongly 
interacting system like superfluid 4 He (for reviews, see Refs. [64} 165]). 

Bose-Einstein condensation manifests itself in Feynman's spacetime ap- 
proach by the formation of so-called cooperative exchange rings, where in- 
dividual particle trajectories hook up to form longer loopsP^ An exam- 
ple is given in Fig. [TH Since Xpm(?i/3) — x 2 (0), xp( 2 )(/i/3) = x 3 (0), and 
Xp(3)(fi/3) = Xi(0), the three particle trajectories now form a single loop, 
wrapping the time axis three times. Generally, a particle in a composite 
exchange ring moves in imaginary time along a trajectory that does not 
end at its own starting position, but at that of another particle. Whence, 
although the initial and final configurations are identical, the particles in 
a composite ring are cyclically permuted. They thereby loose their identity 
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P(3)P(1)P(2) 




1 2 3 



x 

Fig. 16. Worldlines of three nonrelativistic particles on a hypercubic lattice. The vertical 
direction stands for the imaginary time dimension, while the horizontal direction stands 
for all the space dimensions. The particles move synchronously from the bottom to 
the top of the lattice, executing a Brownian random walk in space. Periodic boundary 
conditions are chosen in the time direction. After traversing the lattice, the particles are 
cyclically permuted 1 — > 2 — > 3 — > 1. 



and become indistinguishable. This is the essence of Bose-Einstein conden- 
sation, which at finite temperature can occur in space dimensions larger 
than two. 

The connections between geometrical properties and critical expo- 
nents spelle d out in these notes apply equally well to Bose-Einstein 
condensation!^^ Specifically, the logarithm of the partition function 
P44p takes close to the condensation temperature again the form ([57)l. 
with the worldline loop distribution £b now denoting the number density of 
worldlines wrapping around the time axis b times. In the high-temperature 
phase, the loops have a finite line tension, so that only small loops winding 
around the imaginary time axis only once or at most a few times exist. Upon 
approaching the condensation temperature from above, the line tension van- 
ishes and loops with arbitrary large winding numbers appear { S6 | 27 | s mce 
the winding number of a given loop corresponds to the number of cyclically 
permuted particles in an exchange ring in Feynman's theory, arbitrary large 
rings appear. The particles contained in large rings (long loops) are part of 
the condensate. 

For a free Bose gas in 2 < d < 4 space dimensions, $ = and the fractal 
dimension of the particle trajectories is that of a Brownian random walk, 
D = 2, so that r = d/2 + 1, while a = (d - 2)/2p3 Being related to the 
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thermal critical exponents as before (see Sec. 13. 3p . these exponents yield 

a = d ^ P = h ^ = ih> v = ^ ? ' = ' (148) 

which are precisely the critical exponents of the spherical model in d space 
dimensions. This model corresponds to taking the limit N — > oo in the 
O(N) spin 

model |67|68| Despite being nonintcracting, a free Bose gas is not 
in the universality class of the Gaussian model. The nontrivial exponents 
(|148[) derive from the constraint that the total number of particles is fixed. 
Without this constraint, which is relevant when considering Bose-Einstein 
conde nsation in a free Bose gas at constant pressure, Gaussian exponents 
In four space dimensions, corresponding to the upper critical 
dimension and above, i.e., d > 4, the exponents of the spherical and Gaus- 
sian models coincide. 

The critical properties of Bose-Einstein condensation in an interacting 
nonrelativistic Bose gas are given by the N = 2, if 4 universality class. To 
state this more concretely, the equilibrium thermal critical properties of 
Bose-Einstein condensation in such an interacting system in three space 
dimensions, say, can be computed from the theory ([3]) with d = 3, now 
considered a classical theory in three space dimensions. The connection 
between the nonrelativistic quantum theory and the classical 0(2) model 
is as follows. 

Consider the square end-to-end vector [x r (r) — x r (0)] of the rth par- 
ticle. For a free Bose gas, with the partition function (|146j) . its average is 
readily calculated with the result in d space dimensions 

([x ? .(r)-x r (0)] 2 \ =d-r. (149) 

This expression shows that the worldlines of free nonrelativistic particles 
can be interpreted as trajectories traced out in space by a random walker 
taking steps of typical size 

(150) 

during each imaginary time interval t = hf3. Here, At is the de Broglie 
thermal wavelength 




(151) 

V m 

The step size of the coarse-grained random walk is small when the ther- 
mal wavelength is small, i.e., at high temperatures {(3 small) and for large 
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particle masses. These conditions correspond to the classical limit. Bose- 
Einstein condensation sets in when the step size becomes on the order of 
the interparticle distance l/n 1 ^. Setting a = \jn x l d leads to an estimate 
of the condensation temperature of a free Bose gas 

2ntf ' - N 2,d 



k » To = ^[wm) ' (152) 

with £ the Riemann zeta function 

c(*) = £ p (153) 

{=1 

and where the particle number density n is assumed to be given. If this 
density is small or if the mass of the particles is large, the condensation 
temperature is low. 

The inclusion of a self-interaction in the free (3+l)-dimensional nonrel- 
ativistic quantum field theory changes the fractal structure of the world- 
lines, as they acquire a steric repulsion now, but leaves the basic picture 
unchanged. Consider a grand canonical nonrelativistic ensemble 

Z = Y±f*Z u (154) 

with \i the chemical potential and Zi given in Eq. (I146|) . The coarse-grained 
random walks of the interacting theory in three-dimensional space can be 
identified with the high-temperature graphs of the three-dimensional 0(2) 
spin model. Each set bond corresponds to a particle wrapping the time axis 
once. On a cubic lattice, the smallest closed graph involves four bonds, so 
that translated back to the nonrelativistic quantum theory, exchange rings 
smaller than four are not visible in the spin model. In the high-temperature 
phase of the spin model, only a few small high-temperature graphs are 
present. This reflects the fact that in the nonrelativistic quantum theory 
most worldlines wrap the imaginary time axis only once and that very 
few larger exchange rings are present at high temperatures. A large single 
closed high-temperature graph of b bonds amounts in the nonrelativistic 
quantum theory to a large exchange ring involving b particles. The prolif- 
eration of high-temperature graphs in the spin model corresponds in the 
nonrelativistic quantum theory to the proliferation of worldlines, wrapping 
arbitrary many times around the imaginary time axis, and signals the onset 
of Bose-Einstein condensation. 



February 6, 2008 22:45 



WSPC/Trim Size: 9in x 6in for Review Volume 



3ding 



50 W. Janke and A. M. J. Schakel 

5.3. Summary 

The formalism developed in these notes is shown to be easily adapted to 
describe charged systems as well as Bose-Einstein condensation. 

6. Dual Theories 

As was first pointed out by Helfrich and Mullerj^ the high-temperature 
graphs of O(N) spin models describe in addition to the particle trajectories 
at the same time a second set of physical lines. In this section, we briefly 
discuss examples in d — 2, 3 and 4, respectively. 

6.1. Peierls Domain Walls 

Consider the Ising, i.e, the 0(1) model on a square lattice. According to 
the famous Kramers- Wannier duality!^ th e hig h-temperature graphs of the 
model form line- like Peierls domain walls}^ separating geometrical spin 
clusters of opposite orientation on the dual lattice. Each bond in a high- 
temperature graph cuts a nearest neighbor pair of anti-parallel spins on the 
dual lattice in two. The link connecting the two anti-parallel spins on the 
dual lattice is perpendicular to the high-temperature bond on the original 
lattice. On an infinite lattice, the Kramers- Wannier duality implies that 
observables calculated at an inverse temperature (3 in the original Ising 
model can be transcribed to those of the dual model at an inverse tem- 
perature p. The duality map interchanges the low-temperature and high- 
temperature phases. The relation between the two temperatures is readily 
established by noting that an occupied high-temperature bond represents 
a factor K — tanh/3. The anti-parallel spin pair on the dual lattice, that is 
cut by the high-temperature bond, carries a Boltzmann weight exp(— 2/3), 
so thatPS 

tanh/3 = e~ 2/§ , (155) 

or sinh2/3 = 1/ sinh2/3. The critical temperature j3 c follows from this rela- 
tion by setting (3 = (3, with the result 

/? c = ln(l + V2)/2 = 0.44068679 • • • . (156) 

By the Kramers- Wannier duality, the fractal dimension of the line- 
like Peierls domain walls at criticality coincides with the one of the high- 
temperature graphs, i.e., D = 
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6.2. Vortex Lines 

We next consider the three-dimensional 0(2) model. The fractal dimen- 
sion of the corresponding high-temperature graphs has been estimated in 
a recent Monte Carlo stud y^ as D — 1.7655(20) at the critical point. 
By a duality mapP^ ^^ these high-temperature graphs describe at the 
same time the magnetic vortex lines of the three-dimensional Higgs model 
(|132p . Vortices constitute the topological defects of the Higgs model and 
are line-like in three dimensions. The duality map again interchanges the 
low-temperature and high-temperature phases of the two models. Given the 
values^ v = 0.0380(4) and v = 0.67155(27) of the 0(2) model obtained 
through Monte Carlo simulations, the relations (|10T|) and (| 1 03|) lead to the 
values 

r = 6.2965(60), = 0.1965(20), a = 0.8434(13), (157) 

for the configurational entropy exponents of closed and open magnetic vor- 
tices, t and respectively, and for the exponent cr, characterizing the 
vanishing of the line tension when the critical point is approached. The 
estimate for ■& was first given in Ref. [38] . 

From the perspective of vortices, the phase diagram of the three- 
dimensional Higgs model is as follows. For convenience, we take the model 
as effective theory of ordinary type-II superconductors in three space di- 
mensions. Whereas in the superconducting phase only a few small vortex 
loops are present, the density and size of the loops increase when the tem- 
perature rises. Loops of all sizes appear at the critical point, and the vacuum 
becomes filled with vortex loops. The disordering effect of the proliferating 
vortices destroys the superconducting phase and drives the phase transi- 
tion to the normal phase. Because of the Meissner effect, the magnetic 
interaction is screened in the superconducting phase and vortex line ele- 
ments experience only a short range interaction which is well approximated 
by a steric repulsion. This is precisely the type of interaction experienced 
by the high-temperature graphs of the 0(2) model. 

The three-dimensional 0(2) model also possesses vortex lines as topolog- 
ical solutions. In contrast to those in a superconductor, the vortices in a su- 
perfluid experience a long range interaction, mediated by the gapless Gold- 
stone mode associated with the spontaneous breaking of the 0(2) symmetry 
in the \ip\ theory. The model whose high-temperature graphs descr ibe these 
vortices is the Higgs model, as follows again from the duality mapP^ ^^ 
The long range interaction experienced by the high-temperature graphs of 
the Higgs model is the Biot-Savart interaction discussed in Sec. 15.11 The 
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duality between the 0(2) and the Higgs models implies that the particle 
trajectories of one model appear as vortex lines in the otherff^ 

It was first suggested by Onsager-^2l that proliferating vortex loops drive 
the superfluid phase transition in liquid 4 He. He envisaged that as the 
critical temperature is approached from below, the vortex loops prolif- 
eratelooplproliferation and thereby disorder the superfluid state, causing 
the system to revert to the normal state. The precise nature of the three- 
dimensional superfluid phase transition was recently inv estig ated from this 
perspective in a recent high-precision Monte Carlo study!^ One of the ob- 
servables considered is the total vortex line density v. By means of standard 
finite-size scaling analysis of the corresponding susceptibility 



the inverse critical temperature (3 C was estimated and shown to be consis- 
tent with the estimate of a previous study directly in terms of the field 
variables!^ However, when percolation observables were considered, such 
as the probability for the presence of a vortex loop spanning the lattice, 
slight but statistically significant deviations from j3 c were found. For all 
observables considered, the vortex proliferation threshold /3 P is larger than 
(3 C . That is, from these observables one would conclude that the vortices 
proliferate too early at a temperature below the critical one. From the du- 
ality map alluded to above, the vortex proliferation threshold is expected to 
coincide with the critical point. It is not clear at the moment, whether this 
discrepancy is physical, or an artifact of the percolation observables used. 
With the percolation threshold taken as an adjustable parameter, reason- 
able estimates were obtained from percolation observables for the critical 
exponents v and (3, consistent with those of the XY model!^ 

Even in systems without a thermodynamic phase transition, the no- 
tion of vortex proliferation can be useful in understanding the phase struc- 
ture. An example is provided by the three-dimensional Abelian Higgs lattice 
model with compact gauge field (see below) In addition to vortices, the 
compact model also features magnetic monopoles as topological defects, 
which are point-like in three dimensions. In the presence of monopoles, vor- 
tex lines no longer need to be closed as they can originate at a monopole 
and terminate at an antimonopole. When the vortex line tension is finite, 
monopoles and antimonopoles are tightly bound in pairs. The part of phase 
diagram where this is the case is called the Higgs region, which corresponds 
to the superconducting phase in the noncompact theory. When the vortex 
line tension vanishes, monopoles and antimonopoles are no longer bound 



X = L 3 ((v 2 )-(v) 2 ), 



(158) 
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in pairs, and instead form a plasma. In the part of the phase diagram 
where this is the case, the system exhibits charge confinement. It is well 
established that one can move from the Higgs region into the confining 
region without encountering thermodynamic singularities The suscepti- 
bility data for various observables define, however, a precisely located phase 
boundary. More specifically, for sufficiently large lattices, the maxima of the 
susceptibilities at the phase boundary do not show any finite-size scaling, 
and the susceptibility data obtained on different lattice sizes collapse onto 
single curves without rescaling, indicating that the infinite- volume limit is 
reached. In Ref. [76] . it was argued that this phase boundary marks the 
vortex proliferation threshold. A well-defined and precisely located phase 
boundary across which geometrical objects proliferate, yet thermodynamic 
quantities remain nonsingular has become known as a Kertesz line. Such a 
line was first introduced in the context of the Ising model in the presence 
of an applied magnetic 

6.3. Monopole Loops 

As a last example, we consider the four-dimensio nal pure compact U(l) 
lattice gauge theory described by the Wilson actional 

S 9 = P [l-cos^]. (159) 

X,fJ,<V 

Here, j3 is the inverse gauge coupling, the sum extends over all lattice sites 
x and lattice directions fi, and 9 X ^ V denotes the plaquette variable 

@x,p.v = k-iiQx.v — (160) 

with the lattice difference operator A V 6 X ^ = 6 x + a i> tfJ ,—9 Xtfl and the compact 
variable Q x ,y, € [ — "")"'); living on the link connecting the lattice site x with 
x + aft. The link variable is related to the continuum gauge field A^{x) via 

= eaA^x). (161) 

The lattice action (| 1 59|) reduces to the ordinary Maxwell action in the 
continuum limit a — > 0, provided one sets 

e = *h- (162) 

As for the compact Higgs model, the pure compact gauge theory possesses 
point-like monopoles in three dimensions. In four dimensions, these defects 
become line-like!^ The monopole loops experience in addition to a steric 
repulsion also a long range Biot-Savart interaction. The field theory in which 
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these monopole lines appear as particle trajectories (high-tem perature 
graphs) is the four-dimensional noncompact Higgs model {T52)| P^ ElZIIa 
That is, the pure compact U(l) lattice gauge theory is dual to the non- 
compact Higgs model. These models possess two phases. The pure compact 
U(l) theory has a confining phase at strong coupling (small 0), correspond- 
ing to the superconducting, phase of the Higgs model., and a Coulomb 
phase at weak coupling (large 0) characterized by a massless photon. Monte 
Carlo simulations 2 ^ show that in the Coulomb phase only a few small 
monopole loops are present. With increasing coupling constant (decreas- 
ing 0), the density and size of the loops increase. At the critical point very 
close to 0c. = 1-0, the monopole loops proliferate and thereby disorder the 
systemPZESl The Coulomb phase gives way to the strong-coupling confin- 
ing phase, where the vacuum is filled with a spaghetti of tangled monopole 
loops. It is worth emphasizing that the monopole loops drive the phase tran- 
sition, i.e. , the m onopole loops proliferate right at the confinement phase 
transitionP2ESl l n showing this, percolation observables of the type dis- 
cussed here have been usedP3 The order of the phase transition is not an 
established fact. In case it is continuous, the relations between the fractal 
structure of loops and the critical exponents laid out in these notes also 
apply here. 

6.4. Summary 

Besides picturing the worldlines of the particles described by the field theory 
under consideration, the high-temperature graphs can have a second inter- 
pretation as topological line defects, such as Peierls domains walls (d = 2), 
vortices (d = 3), or monopole lines (d — 4). The field theory in which such 
line-like configurations appear as topological defects is said to be dual to 
the original one. Both field theories describe the same physical lines, once 
interpreted as particle trajectories, once as line defects. 
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